Astron. Nachr. / AN 333, No. 5, 433-437 (2012) / DOI 10.1002/asna.20121 1682 



Cosmic Flows surveys and CLUES simulations 



<N 

o 

<N 



o 

u 

6 

h 



> 
<N 

in 
o 

<N 



H.M. Courtois 1 2 >* and R.B. Tully 2 

1 University of Lyon; UCB Lyon 1/CNRS/IN2P3; IPNL, 4 rue Enrico Fermi, 69622 Villeurbanne, France 

2 Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA 

Received 2012 Jan 17, accepted 2012 Mar 12 
Published online 2012 Jun 1 

Key words Cosmology : observations, distance scale; Galaxies : distances and redshifts, Local Group 

Cosmic Flows is a program to determine galaxy distances for 30,000 galaxies with systematic errors below 2%, almost 
ten times the number currently known and a five-fold improvement in systematics. The resultant velocity field will pro- 
vide input for constrained local universe simulations: CLUES (www.clues-project.org). The observed and the simulated 
universe are then comparatively studied. This synergy of observations and theory distinguishes the program, and should 
lead to fundamental discoveries regarding the sources of deviations from the expansion of the universe. Specifically, the 
program should give a definitive answer to one of the most outstanding unsolved problem in cosmology: the cause of the 
motion of 630 km/s of our Galaxy manifested in the microwave background dipole. This paper presents current results 
with particular emphasis on the "great attractor" reconstruction. 
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1 Introduction 

The universe has assembled itself into a complex web over 
the course of 14 billion years. Structure has formed due 
to gravitational collapse driven by dark matter. It appears 
that today repulsive dark energy has gained ascendancy over 
dark matter, leading to the termination of collapse on large 
scales. Yet we have little understanding of the nature of dark 
matter or dark energy, the two dominant constituents of the 
universe. The goal of this program is to contribute to an im- 
proved and accurate understanding of dark matter and dark 
energy. 

On scales of groups and clusters, and of their constituent 
galaxies, stars, and planets, collapse has occurred and a mem- 
ory of initial conditions has largely been erased by dissipa- 
tive processes. However, on scales larger than groups the 
primordial contest between cosmic expansion and gravita- 
tional collapse is still being played out, modulated by a 
reappearance of expansion in the vacuum. Dark matter man- 
ifests itself through deviations from the universal expansion. 
We know with precision that deviations can be important: 
the dipole temperature variation in the Cosmic Microwave 
Background informs us that our Galaxy has a deviant mo- 
tion of 630 km/s in a well defined direction (Fixsen et al. 
1996, ApJ, 473, 576). The preferred assumption is that we 
are being attracted by a lot of dark matter, but we do not 
know by how much, and how far away it is. The cause of 
our deviant motion is an outstanding problem in cosmol- 
ogy. There can be no pretence that we understand the com- 
position of the universe if we cannot plausibly explain our 
peculiar motion. 
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The distribution of dark matter can be inferred if we 
have a detailed map of the departures from cosmic expan- 
sion. The radial component of the 'peculiar velocity of an 
individual galaxy, V pec , can be determined if we know a 
good distance, d, for the galaxy since the (easily established)! 
observed velocity of the galaxy, Vobs, is the sum of the pe- 
culiar velocity and a cosmic expansion component, Hq * d 
where Hq is the Hubble Constant: V Q ba = Hq * d + V pec . 
The Hubble Constant can be established by averaging over 
a sufficiently large and distributed sample with distances for 
which < Vpec >= 0. 

There was community enthusiasm for peculiar velocity 
studies during the 1990s reaching a crest with the Cosmic 
I Flows workshop in Victoria, Canada (Courteau et al. eds. 
2000, Cosmic Flows Workshop). There has been a hiatus 
since because it came to be appreciated that a lot more and 
better data were needed. The case will be made here that 
vastly superior data are at hand and, moreover, that more 
sophisticated analytic tools are available. We intend to be 
on the vanguard of a resurgence of activity in cosmic flows 
studies. 

The main observational difficulty is with the measure- 
ment of galaxy distances. Say a distance is measured with 
an rms uncertainty of 15%. For a galaxy with V Q bs = 4,000 
km/s the uncertainty of the cosmic expansion component 
(hence the uncertainty in V pec ) will be 600 km/s, compa- 
rable to our bulk motion with respect to the Cosmic Mi- 
crowave Background. Nearby we do well. For a galaxy at 
600 km/s the error is only 90 km/s. At 10,000 km/s the error 
is 1,500 km/s. Peculiar velocities of individual galaxies can 
be mapped in great detail nearby but at larger distances we 
only get statistical information. The numerical orbit recon- 
struction and simulations that will be discussed play a Clit- 
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ical intermediary role in the interpretation of line-of- sight 
peculiar velocities in terms of 3D motions. Maximum like- 
lihood models can be fit that tailor to low errors nearby and 
in high density regions and high errors at a distance and 
where data are sparse. 

Nowadays there are several good ways to measure galaxy| 
distances but there are limitations to each method. We will 
incorporate all methods that are considered appropriate al- 
though our observational efforts will focus on the only method] 
that can provide a dense grid of measurements over a large 
local volume even in low density regions. Methodologies 
we will incorporate include those based on the Cepheid Period-| 
Luminosity Relation (Freedman et al. 2001, ApJ, 553, 47), 
the Tip of the Red Giant Branch (Rizzi et al. 2007, ApJ, 661, 
815), Surface Brightness Fluctuations (Tonry et al. 2001, 
ApJ, 546, 681), the Fundamental Plane (Colless et al. 2001, 
MNRAS, 321, 277), and Type la Supernovae (Tonry et al. 
2003, ApJ, 594, 1). These procedures can give distances 
to individual galaxies or groups with accuracies of 5-10%. 
However, either they can only be applied to very nearby 
galaxies, or galaxies of early morphological type, or serendip-| 
itously. Our observational campaign is directed at applica- 
tion of the Tully-Fisher Relation (Tully & Fisher 1977, A&AJ 
54, 661), the correlation between the luminosity of a galaxy 
and its rotation rate. This methodology gives individual dis- 
tance measures accurate to 15%, less precise than others but 
it is applicable to most spiral galaxies, the most common 
and widely distributed kind of luminous galaxy. Our pro- 
gram will measure accurate distances to 30,000 galaxies: an 
order of magnitude increase over the present state-of-the- 
art. 

A measurement of a distance with the Tully-Fisher Rela- 
tion requires two observations: one to provide the luminos- 
ity of old stars, inferred to be correlated with the total mass, 
and one to determine the rotation of the galaxy, a kinematic 
response to the total mass. Distance is given by the offset be- 
tween the intrinsic luminosity gauged from the rotation rate 
and the observed luminosity. While the methodology has a 
long history, several game-changing advances are about to 
be realized affecting both measurables. 

2 Radioastronomy surveys 

Galaxy rotation can most accurately be determined by ob- 
serving the HI radio 21cm line. The important new and an- 
ticipated capabilities are provided by the upgraded Arecibo 
Telescope in Puerto Rico, the 100 m Green Bank Telescope 
in the USA, ASKAP, the Australian Square Kilometre Ar- 
ray Pathfinder and, ultimately, the wide field-of-view up- 
grade to the Westerbork Synthesis Radio Telescope, Apertif , 
in the Netherlands. We have been engaged in observations 
with the Arecibo and Parkes telescopes and, most impor- 
tantly, with a Large Program on the Green Bank Telescope 
(http://ifa.hawaii.edu/cosmicflows). As a result we have a 
well established data analysis pipeline and, with the inclu- 
sion of archival material, have assembled a database of 16,000| 
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Fig.l The HI profile for PGC31995 as previously ob- 
served at Nancay (top) and as now obtained with the Cosmic 
Flows Large Program on the Green Bank Telescope (bot- 
tom). 



HI profiles (Courtois et al. 2009, AJ, 138, 1938); see http://edd.ifa.hawaii.edul 
and select the 'All Digital HI catalog. The next giant step 
will be with ASKAP with observations beginning in 2011. 
It is anticipated to detect 200,000 extragalactic HI sources 
within z = 0.23 over three-quarters of the sky. We will use 
our innovative newly established procedures to determine 
rotation rates for all adequately detected sources (8 sigma 
detections of the peaks). The final quarter of the sky will be 
observed with Apertif (Westerbrok Netherlands), providing 
galaxy rotation information for 30,000 galaxies, with high 
density coverage within 10,000 km/s and some information 
as deep as 20,000 km/s. 

3 Surface photometry surveys 

In the past, the photometry necessary to complement the 
kinematic measurements have been carried out at optical 
and near infrared bands by a large number of players using 



© 2012 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2012) 



435 



many telescope, detector, and filter combinations in condi- 
tions with uncertain controls. The results have been unsat- 
isfactory. A simple thought experiment illustrates the prob- 
lems experienced with photometry in the past. Suppose there| 
is a systematic inconsistency of 10% between north and 
south or spring and fall or one literature source and another. 
Such a plausible inconsistency would act on targets at 4,000 
km/s to generate a false flow of 400 km/s. Our expectation 
is that the new optical and infrared surveys will enable sys- 
tematics to be held below 2%. 

The way forward has been demonstrated by two ma- 
jor projects. The Two Micron All Sky Survey (2MASS) re- 
sulted in a catalog of all galaxies brighter than Ks = 13.5 
(Jarrett et al. 2000, AJ, 119, 2498). The Sloan Digital Sky 
Survey (SDSS) observed the northern polar cap in 4 optical 
bands and likewise resulted in a catalog of galaxies to faint 
magnitudes (http://www.sdss.org). In both cases, the photo- 
metric integrity for point sources is outstanding. However 
in the case of 2MASS the exposure times were short and 
the sky background at 2 microns is intense with the con- 
sequence that only the highest surface brightness compo- 
nents of galaxies are registered. In the case of SDSS, the 
pipeline processing had poor control over sky levels with 
large galaxies. It is possible to reprocess and recover ade- 
quate photometry for big galaxies on a case-by-case basis 
but, still, SDSS only covers a quarter of the sky. 

We are involved with two alternative approaches to the 
photometry problem either one providing an adequate solu- 
tion. The first approach derives from two ground-based opti- 
cal surveys, PanSTARRS in Hawaii (http://pan-starrs.ifa.hawaii 
and SkyMapper in Australia (http://msowww.anu.edu.au/skymc 
which between them are covering the entire sky in g,r,i,z 
bands. Both surveys are just beginning and over the course 
of 5 years will generate increasingly deep images of the en- 
tire sky with rigorous photometric quality control. The sec- 
ond solution involves imaging at 3.6 microns from space 
with either of two satellites: Spitzer and WISE. In the case 
of Spitzer, observations are made in pointing mode and al- 
ready 2000 galaxies have been targeted and observations 
have been granted (200 hrs) from 2011. We are also mem- 
bers of the consortium S4G (Kartik Seth) who have already 
observed 3000 TF galaxies. On the other hand WISE, the 
Wide-Field Infrared Survey Explorer, has now completed 
imaging of the entire sky. Sky background which plagues 
ground observations in the infrared is not a problem from 
space. Spitzer data are extremely deep and robust for our 
purposes. WISE data are less deep but adequate. In the case 
of WISE, there is not yet an extended source processing 
pipeline so there is work currently done in this direction. 

4 Databases 

The material from the large programs that have been de- 
scribed and more need to be managed in ways that go be- 
yond the requirements of the individual large radio, infrared, 
and optical surveys described above. We have experience 
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Fig. 2 Photometry results from University of Hawaii 2 . 2m 
observations of the galaxy PGC31995 (whose HI profile is 
illustrated in FigQ]). Anti-clockwise from lower left: image 
of the galaxy with stars masked and isophotal contours su- 
perimposed; blow-up of center; surface brightness as a func- 
tion of radius; two panels on different scales of ellipticity 
and position angle. 



with the co-management of an important extragalactic database| 
at the University of Lyon : HyperLeda (http://leda.univ-lyonl .fr).| 
The new database we are developing (Tully et al. 2009, 
AJ, 138, 323) is called EDD, the Extragalactic Distance 
Database (http://edd.ifa.hawaii.edu). It serves three primary 
functions. First, it provides a repository for data that we 
have obtained through observations or extracted from archi ves| 
and reprocessed. The most consequential current asset in 
EDD is the 'All Digital HI catalog of 16,000 consistently 
processed HI profiles for 14,000 galaxies (Courtois et al. 
2009, AJ, 138, 1938). 
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The second primary function of EDD is to assemble the 
many distance scale building blocks so they can be eval- 
uated and linked. The database consists of catalogs, some 
extracted from the literature and some constructed by us, 
but in every case the line entries pertain to a specific galaxy 
identified by its 7 digit PGC name (Paturel et al. 1996, ISBN 
2-908288-08-7). Presently there are information on 50,000 
galaxies. EDD is open without password. 

The third primary function of EDD is to facilitate the 
generation of three-dimensional spatial and velocity field 
maps of the nearby universe. The observed building blocks 
must be converted into distances, alternative distances by 
different methodologies must be evaluated for consistency, 
alternative distances to the same galaxy must be reconciled, 
the grouping properties of galaxies is to be understood and 
mean distances to groups are to be found. Then a value of 
the Hubble Constant consistent with the distance scaling 
should be determined, hopefully derived by averaging over 
a large enough volume to be representative of the global 
value. With all this information, line-of-sight peculiar ve- 
locities can be separated from observed velocities. 

5 Constrained cosmological simulations 

Some of our initial experiences have been with two large 



Alternative approaches are explored, they are based on 
the reconstruction of galaxy orbits to todays observed po- 
sitions and velocities from initial conditions that globally 
minimize the action. One procedure developed by colleagues| 
in France that goes by the acronym MAK, Monge-Ampere- 
Kantorovich (Lavaux et al. 2010, ApJ, 709, 483), finds straight-| 
line orbits for very large samples. This method can be used 
to roughly determine the initial conditions of redshift sur- 
veys of hundreds of thousands of galaxies. A second pro- 
cedure called NAM, Numerical Action Methods (Peebles 
et al. 2001, ApJ, 554, 104), is restricted in use to smaller 
samples of order 1000 objects but is able to find physically 
meaningful orbits in the non-linear regimes of high over- 
density with respect to the universal mean. These several 
methodologies are complementary. MAK is useful to de- 
velop a coarse picture on a large scale. NAM is useful on 
small scales where detailed information is available. It can 
be embedded in a MAK reconstruction. The combination 
can constrain the initial density perturbation conditions lo- 
cally. The Wiener filter approach of CLUES provides an al- 
ternative route to a description of the initial conditions. A 
very attractive feature of the Wiener filter approach is that a 
density map can be constructed by extrapolation to of order 
twice the dimensions (ten times the volume) of the domain 
of distance measurements because of coherence in the prop- 
erties of flows on large scales. The initial conditions that can 



scale simulation projects LEGO (http://graal.ens-lyon.fr/LEGQi recovered can serve as the starting point of a constrained 
and 'HORIZON (http://project-Horizon.fr). Progress were N _body and hydrodynamic simulation, the results of which 
made by the development of the skeleton technique to demon-| are t0 be compared with the observed world, 
strate large scale connectivity (Sousbie et al. 2008, ApJ, 
672, LI), the development of procedures for the generation 
of realistic mock catalogs (Sousbie et al. 2008, ApJ, 678, 
569), and the development of a web interface to launch large 
cosmological simulations on a grid of computers in France 
(Depardon et al. 2010, AIPC, 1241, 816). 

Most recently, we have joined the larger international 
collaboration CLUES, Constrained Local UniversE Simula- 
tions (http://www.clues-project.org) with four Pis: S. Gott- 
loeber in Germany, Y. Hoffman in Israel, A. Klypin in USA, 



6 Current results 

For our first generation of results from the "cosmic flows" 
program, we used our currently available catalog of distance 
measurements containing 1,797 galaxies within 3000 km/s: 
Cosmicflows- 1 . The Wiener Filter method was used to re- 
cover the full 3D peculiar velocity field from the observed 
map of radial velocities and to recover the underlying lin- 



andG. Yepes in Spain (Gottloeberet al. 2010, arXiv: 1005.2687^ density field (Courtois, Hoffman et al. 2012 ApJ 744, 



The collaboration is generating a map of structure within 
10,000 km/s where the input is the observed peculiar veloc- 
ity field, the response to the distribution of matter. Peculiar 
velocities and density perturbations are directly related in 
linear theory. Wiener filtering (Hoffman 2009, LNP, 665, 
565) of the observed velocity distribution recovers the 3- 
dimensional density distribution. Observers in the collabo- 
ration provide the observational input to the models. The 
task not only involves providing the velocity field informa- 
tion described above. It also requires making the translation 
between what observers see and modelers create. Observers 
see sparsely distributed galaxies but only the most luminous 
ones at larger distances and not those in regions of obscu- 
ration. Modelers create dense webs of particles unencum- 
bered by losses with distance or a zone of obscuration but 
the interpretation of those particles in terms of galaxies is 
an art. 



43). The velocity field within a data zone of 3000 km/s 
was decomposed into a local component that is generated 
within the data zone and a tidal one that is generated by the 
mass distribution outside that zone. The tidal component is 
characterized by a coherent flow toward the Norma-Hydra- 
Centaurus (Great Attractor) region while the local compo- 
nent is dominated by a flow toward the Virgo Cluster and 
away from the Local Void. A detailed analysis shows that 
the local flow is predominantly governed by the Local Void 
and the Virgo Cluster plays a lesser role. 

The motion identified as 'Local' is in excellent agree- 
ment in both amplitude and direction with the local motion 
identified by Tully et al. 2008 (a motion of the Local Sheet 
of 323 ± 25 km/s toward l gal = 220 ± 7°, b ga i = +32 ± 6° 
or sgl = 80°, sgb = —52°). The dominance of the Local 
Void marginalizes the dynamical role of the Virgo cluster 
on the local scene. The motion of the Local Group with re- 



© 2012 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2012) 



437 




Fig. 3 3D reconstruction using the Wiener Filter method. 
The great attractor infall is seen on the center left of the 
image. 

spect to the CMB is recovered by the present reconstruction 
to within 5% in amplitude and 5 degrees in direction. 

Decomposing the velocity field into local and tidal com- 
ponents, with respect to the data zone sphere of R=3000 
Mpc.h~ 1, we find that the local structure contributes a com-| 
ponent of 304 km/s leaving a residual tidal component of 
382 km/s, both deviating from the CMB dipole by 30° but 
in almost orthogonal directions. 

Additionally a method is being developped to recon- 
struct the cosmological displacement field using the pecu- 
liar velocity of dark matter haloes at z = as a proxy. We 
apply the Zel'dovich approximation in reverse to link the ve- 
locity and displacement fields and to trace the haloes back 
to their initial position. We quantify the error of this pro- 
cedure depending on different properties of the haloes and 
their environment (Doumler et al. 2012 in prep.). We then 
investigate how this method can be applied to realistic radial 
peculiar velocity data. 



ApJ, 404, 441) while velocity field maps provide informa- 
tion that is fine-grained locally and degrades with distance. 

On the observational side the new capabilities for spec- 
troscopy in the HI line with the new radio telescopes and 
photometry with all-sky surveys from the ground at opti- 
cal wavelengths and from satellites in the infrared are trans- 
formative. We can now envisage a sample for distances of 
30,000 galaxies with integrity in the photometric calibration 
at the level of 2%. This is a very high gain over the existing 
situation (of order 4,000 galaxies with some sort of distance 
measured and with systematics at the level of 10%). 

Through a convergence of observational opportunities 
and computational and display capabilities it is now pos- 
sible to make serious progress with one of the outstand- 
ing problems in cosmology: the nature of deviant galaxy 
motions presumably a response to the properties of the 
dark sector. The deviant motions are seen over a wide range 
of scales. Locally, random motions are seen to be remark- 
ably small but there can be abrupt discontinuities. On large 
scales, coherent flows of impressively large amplitude are 
seen. There is still no agreement on the scales and ampli- 
tudes. There are hints that these flows are larger than can 
be understood by the standard Lambda Cold Dark Matter 
paradigm. Researchers and the public are puzzled by the 
legitimate question: Where are the Dark Energy and Dark 
Matter that dominate the cosmos? 
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1 as initial conditions will be the second generation of the 
CLUES constrained simulations. 



7 Conclusions 



Peculiar velocities are signatures of the distribution of dark 
and luminous matter. Collapse timing is a signature of dark 
energy. Velocity field studies are complementary to weak 
lensing, the only alternative method for directly measuring 
the amplitude of dark matter fluctuations on scales that have 
not yet collapsed. Weak lensing provides statistical informa- 
tion on relatively distant structures (Kaiser & Squires 1993, 
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